clear all
clc

s = tf('s')

sys = 1/(s^2+1.6*s+1)

step(sys,'r')
hold on
title('');
xlabel('Tid [s] ');
ylabel('Position [m]');
axis([0 10 0 1.2])

hline = refline([0 0.98])
set(hline,'Color','b')

hline = refline([0 1.02])
set(hline,'Color','b')

line(1,1)